Mon. Not. R. Astron. Soc. 000, [TJEI (2011) Printed 18 October 2012 (MN LM^X style file v2.2) 



Dead Zones In Circumplanetary Discs as Formation Sites 

For Regular Satellites 



csj! Stephen H. Lubow and Rebecca G. Martin 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

o 



o 
O 



w 

43 

Of 
6 

h 



> 

in 
d 

(N 



13 



ABSTRACT 

Regular satellites in the solar system are thought to form within circumplanetary discs. 
We consider a model of a layered circumplanetary disc that consists of a nonturbulent 
midplane layer and and strongly turbulent disc surface layers. The dead zone provides 
a favorable site for satellite formation. It is a quiescent environment that permits 
the growth of solid bodies. Viscous torques within the disc cause it to expand to a 
substantial fraction of its Hill radius (~ 0.4i?n) where tidal torques from the central 
star remove its angular momentum. For certain parameters, the dead zone develops 
into a high density substructure well inside the Hill sphere. The radial extent of the 
dead zone may explain the compactness of the regular satellites orbits for Jupiter and 
Saturn. The disc temperatures can be low enough to be consistent with the high ice 
fractions of Ganymede and Callisto. 
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1 INTRODUCTION 

After a planet forms in a circumstellar disc and be- 
fore its mass reaches a value of Jupiter's mass, 
tidal forces from the planet open a gap in the disc 
jLin fe Papaloizoul 1 1984 Il98rj ; iD'Angelo, Henning fc Kiev! 
2002; iBate et al.l l2003f T The planet continues to accrete 
material from the circumstellar disc through the Rap and 
a circumpla netary disc forms dArtvmqwicz fc Lubowl 



1996; 



1 19991 ; 

D'Angelo, Henning fc Klevl 120031 ; lAvliffe fc Bate! 120091) . 



Lubow. Seibert fc Artvmowicj 



A circumplanetary (spin-out) disc may also form if the 
planet manages to achieve break-up velocities as it contracts 
(|Korvcanskv. Bodenheimer fc Pollack I Il99ll ). The regular 
satellites in the solar system have prograde, nearly circular, 
and nearly coplanar orbits. They are thought to have formed 
in a circumplanery disc, in a nalogy with the formation of 
plane ts in the solar nebula IjPollack. Lunine fc Tittemorel 
Il99ll ). Understanding the structure of these discs is essential 
for explaining regular satellite formation. 

There are several key constraints that various models 
of gaseous circumplanetary discs have attempted to satisfy. 
The regular satellites around Jupiter and Saturn lie within a 
radius of less than 0.06i?H of their Hill (tidal) radii Rh- The 
small extent of this region has been explained to be a conse- 
quence of a compact disc structure or substructure. Various 
studies have tied the compactness to the level of angular mo- 
ment um of accreting material entering a planet 's Hill sphere 
(e.g-. lEstrada et alj|2009l ; IWard fc Canurj|2010t ). The idea is 



that if the accreting inflow has sufficiently low angular mo- 
mentum per unit mass, then the satellites formed with a 
circumplanetary disc will occupy a small region within the 
Hill sphere. 

All hydrodynamics simulations of gas flows about a 
Jupiter mass planet have revealed a disc that is consider- 
ably mor e extended than the G allilean satellites. The recent 
study by lAvliffe fc Bate! (|2009l ) shows that the disc extends 
to about 0.35-Rh- For a planet such as Jupiter that opens 
a gap in the circumstellar disc, gas enters the Hill sphere 
preferentially near the Lagrange points and carries signifi- 
cant angular momentum. However, the main issue is that 
the compactness of a fully turbulent disc c annot be main- 
tained due to the action of vis cous torques (|Ward fc Canupl 
l20ld ; lMartin fc Lubowll2011al ). These torques cause a com- 
pact circum- Jovian disc to rapidly expand in only ~ 100 
planet orbital periods for typical disc parameters. 

Although the angular momentum of the inflowing gas 
is uncertain, the radial extent of the disc is lar gely inde- 
pendent of this quantity (|Martin fc LubowllanTal ). The disc 
expands to a radius where the tidal torques from the cen- 
tral star (Sun) remove angular momentum at the rate it 
is supplied by inflowing gas. This radius is estimated as 
about 0.47?h, where ball istic periodic orbits begin to cross 
ijMartin fc Lubo w 2011a). (Some departures from this value 
occur due to pressure effects.) In addition, no compact sub- 
structure is found to lie within such a circumplanetary disc 
jEstrada et al.ll2009l ; lAvliffe fc Bate]|2009l ; iMartin fc Lubowl 
l2011al ). With this simple alpha disc model, the turbulent 
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viscosity is by assumption a smooth function of radius. The 
disc structure then follows that of a standard smooth accre- 
tion disc. 

Previous work such as ICanup fc Wardl (|2002l ) suggested 
that if the inflowing angular momentum is small enough, the 
flow may penetrate deep into the Hill sphere before initially 
collecting into a disc. If the disk is viscous, its gas compo- 
nent will then spread outward to a large radius. But the 
solids initially delivered with the gas may not be effectively 
coupled to the gas if they rapidly accrete and grow once in 
circumplanetary orbit. The solids then accrete near the com- 
pact region where they were initially delivered, rather than 
tracking the outward viscous expansion of the gas disk. 

There are some issues of concern with this picture. The 
inflow may be largely planar and join the circumplanetary 
disk at its outer edge, as indicated by simulations by Ayliffe 
& Bate (2009). The gas and coupled solids would then not 
be able to penetrate deeply within the Hill sphere. Instead, 
they achieve Keplerian orbits near the outer edge of the pre- 
existing disk as they becomes entrained there. It is possible 
that larger size solid bodies that are decoupled from the gas 
could enter the Hill sphere with sufficiently low angular mo- 
mentum to be captured close to the planet. But it is not clear 
that such a low angular momentum occurs. In this paper, 
we suggest an alternative explanation for the compactness 
of the satellites based on the radial extent of a dead zone, a 
region of low disc turbulence. 

Another constraint applies to the disc temperature. 
This snow line plays an important role in the composition of 
forming satellites. Outside of the snow line, the solid mass 
density is much higher because of water ice condensation. 
The sn ow line occurs at a temperature , T snow , that is around 
170 K (|Havashilll98ll ; lLecar et alj|2006f ). Ganymede and Cal- 
listo, the two outer Gallilean satellites contain a substantial 
fraction of ice (~ 50%), while the two inner satellites contain 
much less. The snow line of the disc is then required to have 
been near the orbits of these outer satellites. The partially 
differentiated structure of Callisto suggests that its ice never 
full y melted and that the snow line was always inside its or- 
bit (|Lunine fc Stevensonlll982h . Since the disc temperature 
increases with mass accretion rate, this constraint limits that 
rate. This requirement implies that the disc accretion rates 
be sufficient low, substantially lower than w hat would be ex- 
pected during the T Taur i accretion phase |Canup fc Wardl 
|2002| ; Ifetrada et"ail 2009). The satellite formation epoch is 
therefore expected to occur late in the life of the solar neb- 
ula, when it has lost a substantial fraction of its mass. 

Another issue involved in the formation of both planets 
and satellites in gaseous discs is the permitted level of tur- 
bulence. Estimates of disc turbulence from simulations and 
properties of observed systems suggest that relatively high 
levels of turbulence a ~ 0.01 — 0.3 are expected in fully 
turbulent discs (e.g., iKing. Pringle fc Livid [2007T I. Larger 
solid bodies may form from the gravitational instability of a 
dense layer of dust that settles near the disc mid-plane (e.g., 
iGoldreich fc Wardlll973l ). Even modest levels of turbulence, 
involving a <C 0.01, can prevent solids from setting to a thin 
enough layer near disc mid-plane f o r gravitational instability 
to operate (|Dubrulle et all 1 19951 : [c uzzi fc Weidenschillind 
120061 ) . An alternative model for the growth of so lids involves 
the c oncentration of dust in turbulent eddies (|Cuzzi et al.l 
2008). This model also adopted a weaker level of turbulence. 



Another effect is that turbulence can cause destructive colli- 
sions among larger solid bodies, thereby preventing growt h 
to the desired size (e.g., Ilda. Guillot. fc Morbidelli 1 120081 ). 
These results suggest that lower values of disc turbulence 
may be more favorable for satellite formation than is occurs 
in a fully turbulent disc. 

In a layered disc , the magnetic turbulenc e (MRI; 
[ Balbus fc H awlcv 1991) is not active at all heights |Gammiel 
Il996h . Magnetic turbulence requires a certain level of ioniza- 
tion for the gas to be well enough coupled to the magnetic 
field. In sufficiently cool disc regions, temperatures are too 
low to provide the needed level of the ionization for magnetic 
turbulence to operate. Instead, the ionization is provided by 
external sources of radiation, such as X-rays or cosmic rays. 
If the disc surface density is high enough, this radiation may 
not penetrate deep enough below the disc surfaces to pro- 
vide turbulence at all heights. A so-called dead zone with 
little or no turbulence develops around the disc midplane. 
Such an environment may be favorable for the survival and 
growth of solid bodies. 

In an earlier paper, we explored a layered disc model for 
the evolution of a circumplanetary disc during a stage when 
the disc experiences accretion rates compa rable to those of 
the T Tauri phase (|Lubow fc Martini l2012h . Such discs can 
sometimes undergo outbursts, analogous to the FU Ori out- 
bursts. In any case, such a disc is too hot for the survival of 
icy satellites at the locations of Ganymede and Callisto. 

We explore in this paper the viability of a layered disc 
model as an environment for satellite formation in later 
stages of its evolution. We do not attempt to model the 
formation of satellites in such a disc. Instead, we show that 
in this stage and for certain plausible disc ionization pa- 
rameters, a layered disc model satisfies the requirements de- 
scribed above. In particular, a compact high density disc 
substructure develops in a low turbulence dead zone whose 
extent matches that of the regular satellites. The disc tem- 
perature is low enough for the outermost Galilean satellites 
to lie outside the snow line. 

In Section [5] we describe the results from our circum- 
planetary disc model and how it may explain some features 
of the formation of the Galilean satellites. In Section [3] we 
discuss some implications and limitations of the model and 
in Section [4] we draw our conclusions. 



2 CIRCUMPLANETARY DISC MODEL 

In this Section we first describe our layered circumplanetary 
disc model that is truncated by tides from the star and then 
we discuss the results. 

2.1 Model Description 

W e use the layered disc model initially described 
in lArmitage. Livio. fc Pringld (l200lf) and f u rther d evel- 
oped in IZhu et al.l (I2009T ). iMartin fc Lubowl (1201101 ) and 
iLubow fc Martini l|2012l ). Material in the circumplanetary 
disc orbits the central planet of mass M p at Keplerian an- 
gular speed Q(R) — \jGM v j R? for radius R from the 
planet. Material is added to the disc at rate Mmfaii- The 
disc has a total surface density T.(R,t), midplane temper- 
ature T c (R,t), and surface temperature T c (R,t) at time t. 
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The MRI turbulent surface layer (which we call the active 
layer) has surface density E m (ii, t), temperature T m (R,t), 
and tur bulent viscosity ^ m (T m ,J?, t) that is parametrised 
with the lShakura fc Sunvaevl (1 19731 ) a m parameter. Cosmic 
rays and/or X-rays are assumed to be able to provide suf- 
ficient ionization for MRI to operate to a constant surface 
density depth E cr it- If the total surface density is greater 
than this critical value, E > E cr i t , then E m = E C rit- For 
smaller surface density, E < E cr it, the disc is fully MRI ac- 
tive and E m = E. 

A complementary midplane layer exists if the total sur- 
face density is larger than the critical, E > E cr i t . It has 
surface density E g = E — E cr i t and midplane temperature 
T c . This layer is MRI active if the midplane temperature is 
greater than the critical, T c > Tent, where T cr it is the tem- 
perature for sufficient thermal ionisation. However, for lower 
temperatures, T c < Tcrit, a dead zone forms. In the model, 
there is no turbulence in the dead zone unless it is massive 
enough to be self-gravitating. (In reality, some turbulence is 
expected du e to the vertical propag ation of waves from the 
active layer (|Fleming fc S tone 2003).) The condition for self 
gravity is taken to be that the Toomr e parameter is smaller 
than its critical value, Q < Q C rit = 2 (|Toomreill96"3 ). In this 
second viscosity term is included in the complemen- 
tary layer. 

We consider the planet Jupiter with mass, M p = 1 Mj, 
that is orbiting the Sun, M = IMq, at a distance of a — 
5.2 AU. We solve the accretion disc equations numerically on 
a fixe d mesh that is uniform in log J? with 120 grid point s 
(e.g., lArmitage. Livio. fc Pringlell200ll ; iMartin et aT1l2007h . 
The inner boundary has a zero torque condition at the radius 
of Jupiter, R = 1 Rj, so the mass falls freely on to the planet. 
At the outer boundary we approximate the tidal torque from 
the star with a zero radial velocity boundary condition at 
Rout = 0.4 Rh- Material accretes on to the disc at a constant 
rate Mj n f a n at a radiu s which we take to be 0.33 -Rh (see 
iLubow fc Martini |2012| . for more explanation). Initially, we 
take the surface density and temperature structure of the 
disc to be that of low mass (nonself gravitating) disc with a 
small dead zone. This state represents the disc immediately 
following an outburst. Alternatively, if there is initially no 
mass in the disc, it quickly builds up to the same state. 

The critical surface density that is sufficiently ionised to 
be turbulent, E C rit, is not well determined. When cosmic rays 
are the dominant source of ionisation, for typical parameters 
the surface turbulent layer of a circumstellar disc has been 
estimated to h a ve a surface density of Sent ~ 200 g cm -2 
<|Gammiel Il996l ; iFromang. Terquem fc Balbus! |2002| ). How- 
ever, it is possible that cosmic rays are swept away from the 
disc by means of a turbulent M HD jet or outflow which 
may be present at th e star (e.g. ISkilling fc Stron gj 1 19761 ; 
ICesarskv fc V5Udll978h . In the absence of cosmic rays, X- 
rays from the star may be the dominant source of ioni- 
sation and in this case the active layer is much smaller 
|Matsumura fc Pudritzll2003h . 

More recent works find the inclusion of polycyclic aro- 
matic hydro carbons and dust tends to suppress the instabil- 
ity further (|Bai fc Goodman! 120091; iPerez-Becker fc Chiang! 
l201ll ; lBaill201ll ; lFuiii. Qkuzumi fc Inutsukall201ll ). However, 
in circumstellar discs, such effects produce an accretion 
rate that is too low to account for T Tauri accretion rates 
and this problem remains unsolved. IPerez-Becker fc Chiang! 



|201ll) suggest that in order to explain the T Tauri rates, 
Ecrit > 10 gem -2 , in the circumstellar disc. It is possible 
that circumplanetary discs may have an even smaller ac- 
tive layer surface density than circumstellar discs, since they 
may be shielded from the ionising radiation of the central 
star by the shadowing effects of the inner circumstellar disc. 
Because of these uncert ainties, we follow the work of Ar- 
mitage et al (2001) and IZhu et all (|2009h and regard E cr it 
as a constant free parameter. If the critical surface den- 
sity is large, then the entire disc will be fully turbulent. 
If it is small, then the dead zone will be very extensive. 
We consider a range of active layer surface densities with 
lgcm -2 s: E C rit < 100 gem" 2 . 

Generally, we take the value of the temperature for 
sufficient thermal ionization for MRI to operate to be 
Tcrit = 800 K. At this temperature the ionisation fraction 
increases exponentially with t emperature due to the colli- 
sional ionisation of potassium (jUm cbavashi 1983). The vis- 
cosity Q m parameter associated with MRI is found to be 
> 0.01 from MHD simulations, but depends on numerous 
parameters such as the resolution, stratification, and treat- 
men ts of small scale dissipation and radiation transport 



e.g. Hartmann et al.lll998l ; lFromang et al.l20o"7l ; lGuan et al.l 



120091 ; iDavis et al]|200ST ). Similarly obser vations of FU Ori 
suggest that a m « 0.01 l|Zhu et al.ll2007l ). However, obser- 
vations of X-ray binaries and dwar f novae suggest a m w 
0.1-0.4 jKing. Pringle fc Lividl2007h . Armitage et al (2001) 
take a m = 0.01 and Zhu et al.l (|2010h consider both 0.01 and 
0.1. With the uncertainty in the value, we also consider cases 
with both Q m = 0.01 and 0.1. 

During times of planet formation, the accretion rate 
on to the circumplanetary disc, Mjnfail, is of order 
the overall circumstellar disc accretion rate (iBate et al.l 
120031 ; ILubow fc D'Angelol [20061; lAvliffe fc Bate! 120091) that 



is of order 10~ 8 Mp vr~* l|Valenti. Basri fc Johns! 1 19931 : 
iHartmann et al.lll998l ). typical for the T Tauri phase. During 
this phase, disc temper atures are too high to permit the sur- 
vival of ice in Callisto (|Canup fc W ard 2002]; lEstrada et al.l 
2009). Consequently, we consider the evolution in the later 
stages of disc evolution, when the accretion rate has dropped 
by an order of magnitude or more. We consider accretion 
rates in the range 10~ n — 10~ 9 Mq yr~* near the end of the 
disc lifetime. 

If the circumplanetary disc density in the dead zone 
grows sufficient large during the T Tauri stage, the disc be- 
comes gravitationally unstable. Self-gravity generates tur- 
bulence that raises the disc midplane temperature (e.g., 
lLodato fc Rice! [20041 ) . If the temperature exceeds the crit- 
ical value Tcrit for sufficient thermal ionization for MRI 
to operate, MRI can set in abruptly and lead to a 
gravo-magneto accretion outburst in a circumplanetary disc 
IILubow fc Martinll2012T ). These conditions involve high disc 
temperatures. The dead zone structure is disrupted and 
dead zone mass is accreted on the central planet. The pro- 
cess repeats and can be d escribed as a limit cycle of the 
gravo-magneto instability (jMartin fc Lubowl l2011bl ) . Con- 
sequently, we regard the T Tauri stage as unfavorable for 
satellite formation or survival. 

As the disc disperses near the end of its lifetime, the 
accretion rate Minfaii decreases and outbursts become less 
likely. We are most interested in models that have a long 
outburst interval timescale. The outburst interval timescale 
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increases with decreasing accretion rate. If the mass accre- 
tion rate decreases substantially after an outburst, then the 
outburst may not recur. The lifetime of the circumplanetary 
disc is assumed to be similar to that of the circumstellar disc, 
a few 10 6 yr. A simple estimate for the outburst timescale 
is the timescale for the disc to become self gravitating. This 
timescale can be estimated as the time to build up a mass of 
MpH/R, where H — c s /fi is the disc scale height and c s is 
the sound speed. This mass is at most around a few tenths 
of a Jupiter mass. The outburst interval timescale is of order 



tn 



2.2 Model Results 



A/ p H 



R' 



(1) 



Table Q] summarizes simulation results for a range of pa- 
rameters that cover some plausible values for the viscos- 
ity parameter a m and the critical surface density for ion- 
ization Ed-it. In Fig. [1] we show the surface density and 
temperature distributions for models with a m = 0.01 and 
Mnfaii = 10" 10 M© yr" 1 . Model R8 is fully turbulent and 
this steady solution holds for £ cr it > 27.1 gem -2 . Model R9 



has critical surface density S c 



10 g cm and hence con- 



tains a dead zone. Within the dead zone, the surface density 
increases substantially with time. For Model R9 this region 
extends to radius i?dead = 0.034 _Rh- This location is close to 
the orbit of Callisto. Such a model may explain why Callisto 
has a large mass and angular momentum, yet no satellites 
are found outside of its orbit. Beyond the dead zone, less 
material is available to form satellites and the turbulence 
may disrupt the growth into larger bodies. 

At later times, the higher disc mass in the dead zone 
is in the r ange of the so-called mini mum mass sub-nebula 
(MMSN) (|Lunine fc Stevenson! Il982^ . However, the model 
here is fundamentally different from a static MMSN that 
was initially envisioned by some earlier work, because it con- 
siders a continuously supplied disc. The accretion time scale 
in this paper is regulated by the inflow rat e to the disc, as 
in th e original continuous inflow concept of |Canup fc Ward! 
120021 ). 

The disc structure we obtain is similar to the gas rich 
disc model envisione d bv lMosqueira fc Estrada! (|2003al ) and 
lEstrada et al.l (|2009h . In their model, Callisto lies just out- 
side the dense inner disc, the dead zone in the current 
model. They sug gested that Callisto's inferred lo ng for- 
mation timescale (jStevenson. Harris fc Lunindll986h , which 
caused its partial differentiation, results from the rela- 
tively slow delivery of solids in the outer disc, due to 
its lower density. The model also describes Ganymede as 
being formed more rapidly within the dense disc. These 
different environments provide a possible explanation for 
the Callisto- Ganymede dichotomy that describes the dif- 
ferences between these tw o outermost Gallilean satellites 
(ILunine fc StevensorJ[l982^ . 

The mass of the fully turbulent disc in Model R8 is 
4.3 x 10~ 5 Mj whereas the mass of the disc with a dead 
zone in Model R9 at a time of 10 s yr is 9.1 x 10" 3 Mj. The 
dead zone model has a lower temperature distribution than 
the fully turbulent model. For example, the snow line loca- 
tion changes from 0.032 _Rh in the fully turbulent model to 
0.017 Rh when a dead zone is introduced. The latter snow 



line location is in the range expected for the survival of ice 
in Callisto. The outburst intervals for Model R9 at constant 
accretion rate are also fairly long. For model R9, H/R < 0.2 
and so the outburst timescale is around 2 x 10 6 yr. 

With a higher accretion rate of 10 -9 M© yr -1 , the disc 
builds up mass very quickly and so the outburst time scale 
is fairly short. In addition, it is difficult to simultaneously 
meet the constraints on the snow line location and the com- 
pactness of the dead zone. The models that best meet these 
requirements, Models R9 and R12, have the intermediate 
accretion rate M = 10~ 10 M Q yr" 1 . 

We note that because of the dead zone, the vertically 
averaged a in the disc is much lower than the value in the 
active layer of 0.01 or 0.1. In some ways t his vertical aver- 
age ap proximates the a value adopted by ICanup fc Ward! 
(2002). However, since that model does not have a dead 
zone, there is no accumulation of mass that produces the 
compact disc substructure. Fully turbulent models with 
lower accretion rates of 10" 11 M© yr" 1 (R13 and R14) have 
low eno ugh temperatures for icy satellite formation (as 
found bv lCanup fc Wardll2002l ). 



3 DISCUSSION 

There are several approximations in the layered disc model 
we have considered that should be impr oved on in future 
work. Some of these were dis c ussed in iMartin fc Lubowl 
( |2011bh and iLubow fc Martini (|2012t ). For example, we 
have approximated the surface density in the active layer 
as a constant with radius in the disc. An alternative 
method of finding the dead zone is with a critical mag- 
netic Reynolds number (e.g. Fromang. Terquem fc Balbusl 



120021 : iMatsumura fc Pudritj 120031 '). In this approach, the 
surface density of th e turbulent layer may vary in radius 
l|Martin et al.l l2012al ). However, there are currently still 
problems with this model, as it is difficult to explain T 
Tauri accretion rate s with a low critical surface density 
jMartin et alj|2012bh . 

The infall accretion rates we considered are constant 
in time. However, they are likely to decrease as the flow 
through the circumstellar disc decreases. As the accretion 
rate decreases, the active surface layer occupies an increasing 
fraction of the vertical extent of the disc. When the accretion 
rate drops to values of < 10 _11 M© yr -1 , the disc will become 
fully turbulent with the adopted parameters (Models R13 
and R14). The disc temperatures remain low enough for the 
survival of ice in Callisto. 

The inferred density of the small inner Jovian satel- 
lite Amalthe a suggests that contai ns a substantial fraction 
of water-ice (I Anderson et al.l 12005). Accretion disc models, 
such as the one described here, predict that temperatures 
are too high for the water-ice to survive and requires some 
other explanation. 

Another open issue is satellite survival in such a high 
density disc against the effects of orbit decay due to disc- 
satellite interactions l|Wardlll997l ; ILubow fc Idall201Ch . The 
Type I migration times cale in the dead zone of model R9 
is quite short ~ 10 2 yr l|Canup fc Wardll2002l ). It is possi- 
ble that the satellites instead undergo the potentially slower 
Type II migration, provided that they can open gaps. Gap 
opening in a low viscosity disc is dominated by effects of gas 
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Model 


Mnfall 






.Rdcad 




Mdi sc (t = 10 5 yr) 


Outburst Timescale 




(Moyr- 1 ) 




(gem 2 ) 




(Ru) 


(Mj) 


(yr) 


Rl 


10- 9 


0.01 


> 90.3 


- 


0.094 


1.9 x 10~ 4 


Steady 


K2 


10 -9 


0.01 


40 


0.023 


0.084 


6.7 x 10" 2 


1.4 X 10 5 


R3 


10" 9 


0.01 


10 


0.32 


0.019 


1.0 x 10" 1 


7.8 x 10 5 


R4 




0.01 


\ 


0.33 




1.1 x lO - 1 


8.3 x 10 5 


R5 


10- 9 


0.1 


> 15.6 




0.056 


2.5 x 10- 4 


Steady 


R6 


lO" 9 


0.1 


10 


0.0097 


0.052 




> 10 5 


R7 


10- 9 


0.1 


1 


0.32 


0.011 


1.0 x lO" 1 


> 10 5 


R8 


10 -io 


0.01 


> 27.1 




0.032 


4.3 x 10" 5 


Steady 


R9 


10 -io 


0.01 


10 


0.034 


0.017 


9.1 x 10- 3 


> 10 6 


RIO 


10- 10 


0.01 


1 


0.32 


0.0023 


1.1 x 10~ 2 


> 10 6 


Rll 


10- 10 


0.1 


> 4.7 




0.021 


4.5 x 10" 6 


Steady 


R12 


10- 10 


0.1 


1 


0.035 


O.Oll 


9.8 X 10- 3 


> 10 6 


R13 


lO" 11 


0.01 


> 8.1 




0.011 


7.7 x 10" 6 


Steady 


R14 


10- 11 


0.1 


> 1.2 




0.0092 


7.8 x 10" 7 


Steady 



Table 1. Column 2 contains the infall accretion rate on to the circumplanetary disc, column 3 contains the viscosity a m parameter, 
column 4 contains the critical surface density in the active layer that is ionised by cosmic rays/X-rays, column 5 contains the outer radius 
of the dead zone, -Rdcadi if it exists, column 6 contains the snow line radius and column 7 contains the total mass in the disc at time 
t = 10 5 yr. If the disc has a dead zone, then the mass of the disc increases linearly in time. But if there is no dead zone, then the mass 
in constant. Finally, column 8 contains estimates of the timescales for the gravo-magneto outbursts. 
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Figure 1. Left: The surface density of the disc in model R8 (dashed line) and R9 (solid lines) at times t = 10 2 , 10 3 , 10 4 , 10 5 and 
10 6 yr in order of increasing height in the plot. As the dead zone gains mass, its surface density grows in time. Right: The midplanc 
temperature of the disc in R9 (solid line) and R8 (dashed line). The temperature in model R9 does not change in time even though the 
surface density does. The dotted line shows the snow line temperature, T snow = 170 K. The solid circles in both plots show the radial 
locations of the Galilean satellites. 



shocks caused by the steepening of density waves ijRafikovl 
120021 : lYu et al.ll20lbT ). The gap opening condition in this 
case suggests that Callisto and Ganymede can open gaps 
for Model R9, but only for £ < 10 4 gcm -2 . 

The nature of gap opening and migration in a layered 
disc has not been examined. The applicatio n of the viscous 
gap opening criterion (e.g., equation 5 of ICanup fc Wardl 
120021 ) to the active layer alone suggests that gap opening 



does not occur. On the other hand, its application to a ver- 
tical averaged a suggests that gap opening may occur. There 
are some other possibilities. If Callisto resides just out- 
side the dead zone, it may be subjec t to outward Lindblad 
torque s that stall its migration there (Mosqueira fe Estrada! 



l2003bl : iMatsumura, Pudritz fc Thommesl 120071 ). Trapping 
may also be possib le at the inner edge of the dead zone 
jKretke et al. | [2009). We also note that outward corotation 
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torqu es could affect migration |Paardekooper fc Mellemal 
120061 ). since the dead zone region in Fig. 1 has a negative 
radial entropy gradient that is required for this effect to op- 
erate. On the other hand, coorotation torques can become 
quite weak and saturate in a low vis cosity environment (e.g., 
IWardlll99ll;lQgilvie fc Lubowll2003l ). 

ICanup fc Wardl (|2002l l point out other constraints that 
should be considered. One is that if satellites open gaps in 
the disc, then satellite eccentricities should be excited by 
first order Lindblad resonances. The excitation of eccentric- 
ity of plan ets that open gaps has been an active area of ex - 
ploration (|Goldreich fc Sad 120031 ; lOgilvie fc Lubowl I2003T ) . 
The simulation resul ts appear to be sensitive to how clean 
and large the gap is. IPapaloizou. Nelson, fc Masset I (|200lh 
found no eccentriciy growth for a lMj which produces 
a fairly clean gap, but g r owth at much higher masses. 
iD'Angelo. Lubow fc Bate! (2006) found some eccentricity 
growth occurred lMj, but was quite weak compared to 
growth rates for higher mass planets. In any case, at late 
times when Mi n f a n < 10 _11 Moyr _1 , the disc surface den- 
sity in the satellite region will drop to less than critical values 
E(i?) < E CI .it. The disc is then fully turbulent and gap open- 
ing will not occur. The available gas at this stage may be 
sufficient to damp eccentricities developed in a prior open- 
gap stage. 

The disc structure produced by this model is 
somewhat similar to the g as r ich model de s cribed in 
iMosqueira fc Estradal l|2003af ) and lEstrada et al l <|2009l ). In 
that model, the dense inner disc is the relic of a disc formed 
prior to gap opening. The low density outer disc is formed 
by higher angular momentum gas accreted after gap open- 
ing. It is not clear why the inner disc would not accrete onto 
the planet. We expect that a significant fraction of the mass 
of Jupiter was acquired after gap opening through disc ac- 
cretion. In the present model, accreting gas builds up mass 
in a small inner region due to the presence of the dead zone. 

The dead zone is assumed to be free of turbulence. 
Suppose some independent mechanism produces turbulence 
there that we do not take into account, such as a hydro- 
dynamic instability not involving MRI or self-gravity. The 
dead zone would then achieve a steady state flow that would 
limit the mass growth of the dead zone as seen in Fig. [1] For 
the parameters adopted in Model R9, reaching the highest 
plotted densities requires that ay ~ 1CP 8 . 



4 CONCLUSIONS 

We have described some possible proto-satellite environ- 
ments that occur in discs with dead zones. A dead zone 
will occur if external r adiation is insu fficient to fully ionize 
cool regions of a disc (|Gammidll996h . The dead zone pro- 
vides a quiescent environment for the survival and growth 
of solid bodies into satellites. As the dead zone gains mass, 
it can provide a high density, compact substructure within 
the circumplanetary disc (see Fig. 1). The regular satellites 
of Jupiter and Saturn lie within a small fraction of their Hill 
radii. A fully turbulent simple alpha disc extends to much 
larger radii and has a smoothly varying density structure. 
The compactness of a high density dead zone provides a pos- 
sible explanation for the small radial extent of the regular 
satellites. For accretion rates appropriate to the late stages 



its evolution, the circum- Jovian disc can be cool enough to 
permit the survival of the ice in Callisto. 

In this paper we have only considered possible disc 
structures that can arise in circumplanetary discs with dead 
zones. More work is required to explore the growth, survival, 
and evolution of satellites in such discs. 
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